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ABSTRACT 

We investigate the pulsational stability of massive (M > 120 Mq) main sequence stars of 
a range of metallicities, including primordial. Population III stars. We include a formulation 
of convective damping motivated by numerical simulations of the interaction between con- 
vection and periodic shear flows. We find that convective viscosity is likely strong enough 
to stabilize radial pulsations whenever nuclear-burning (the e-mechanism) is the dominant 
source of driving. This suggests that massive main sequence stars with Z < 2 x lO"-' are 
pulsationally stable and are unlikely to experience pulsation-driven mass loss on the main se- 
quence. These conclusions are, however, sensitive to the form of the convective viscosity and 
highlight the need for further high-resolution simulations of the convection-oscillation inter- 
action. For more metal-rich stars (Z > 2 x lO""*), the dominant pulsational driving arises due 
to the /(--mechanism arising from the iron-bump in opacity and is strong enough to overcome 
convective damping. Our results highlight that even for oscillations with periods a few orders 
of magnitude shorter than the outer convective turnover time, the "frozen-in" approximation 
for the convection-oscillation interaction is inappropriate, and convective damping should be 
taken into account when assessing mode stability. 
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1 INTRODUCTION 

Several rounds of investigations over the course of the last sev- 
enty years have shown that massive main sequence stars above a 
critical mass near 100 Mq are vibrationally unstable in their fun- 
damental radial mode (e.g. lLeciouxlll94lUSchvyarzschild & Harml 



ll959l : ISimon&Stotherslll970l;IZiebarthlll970l) . In these stars, the 
low density contrast between the stellar core and envelope allows 
the fundamental mode to reach large amplitude near the core, where 
it may readily couple to the highly temperature sensitive nuclear 
energy generation (the 6-mechanism). This results in linear insta- 
bility with growth times shorter than the stellar evolutionary time. 
The resulting non-linear evolution of this instability, which several 
authors proposed could lead to significant pulsation-driven mass 
loss, remains an unsolved problem despite consider able effort (e.g. 
I Appenzelleill 1970l: Iziebarthi 1 970l: |Papalo"izoull 1 973h . 

The introduct ion of OPAL opacities in stellar models in 
the early nineties d Rogers & IglesiasI Il992h enhanced the linear 
instability of mass i ve sta rs of approximately solar metallicity. 
iGlatzel et all j 19931 ; 1 19991) found that the enhanced opacity in the 
stellar envelope due to transitions in heavy element isotopes can 
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produce strong radiative (/e-mechanism) and strange-mode driving. 
Today, these mechanisms are generally believed to be the most im- 
portant sources of fundamenta l mode instability in massive solar- 
metallicity stars ( lGlatzelll2005h . 

In the investigations to date, convection has been treated in 
the "frozen-in" approximation, in which the interaction between 
pulsation and convection is ignored entirely. The validity of this 
approximation is unclear and needs to be explicitly quantified. 
The treatment of convection-pulsation interaction as an enhanced 
viscosity has a long history rooted in st udies of tidal Q factors 
and the solar 5-m inute oscillations (e.g. iGoldreich & NicholsonI 
ll977l : IZahnlll989l) . Only recently, however, have these theoreti- 
cal esti mates been c alibrated wit h numerical experiments. In par- 
ticular, iPenev et^ I ( l2009bl ; I2OI 1 ) have shown that for oscillation 
timescales of the same order as the outer convective turnover time, 
the interaction is well represented by an anisotropic viscosity which 
scales linearly with the ratio of the oscillation and outer convective 
turnover times (as argued for by Zahn 1989). On timescales shorter 
than about one-third of the outer turnover ti me, these authors argue 
for a q uadratic scaling in accordance with IGoldreich & NicholsonI 
( I1977I) . 

In this paper we reconsider the linear stability of radial os- 
cillations in massive stars, including primordial Population III 
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stars. For primordial stars, the lack of metals excludes opacity- 
driven instabilities in the envelope (i.e., /c-mechanism and strange 
modes). The driving found previously for Population III stars 
above about 120 Mq relies on nuclear-driving in the convective 
stellar core (e-mechanism), with resulting gr owth times that are 
much longer than for opacity-driven modes dSaraffe et alj|200ll : 
ISonoi & Umed a 2011). It is thus particularly important to check 
the effect of convective damping on these weakly-driven modes. 
We also extend our investigation to higher metallicities for ~ 100 
Mq models to examine the effect of convective damping on the sta- 
bility of massive stars of a range of metallicities. 

We begin by describing the equilibrium stellar models used in 
this work in Section|2] In Section|3] we describe our quasi-adiabatic 
linear stability analysis, highlighting the important contributions to 
mode driving and damping. We then describe our primary results 
in Section|4]and discuss their implications in Section|5] 



2 EQUILIBRIUM MODELS 

We have computed evolutionary sequences for equilibrium stellar 
models with initially primordial composition (i.e.. P opulation III 
stars) using the MESA star stellar evolution code jPaxton et al.l 
We have also computed model sequences for an initial 
mass of 120 Mq, and a range of metallicities from Z = 2 x 10"* 
to solar (which we take to be Z = 0.02). All models are non- 
rotating and have mass loss turned off. The former approximation 
is made for simplicity, while the latter is justified for Population III 
stars because the line-driving mechanism for massive star winds re- 
lies on the presence of metal-line opacity in the ste llar atmosphere 
(e.g., lLamers & CassineIIill999l:ICastor et al.ll975l) . For the higher 
metallicity models, since we are interested only in the stability of 
massive stars on the main sequence, the effects of mass loss are not 
critical. 

We use the standard mixing length prescription of 
iBohm-Vitensd ( Il958h and a mixing length of ~ 1.6 pressure 
scale heights for the convective mixing. To examine the effects 
of mixing on pulsational stability, we have tested schemes for 
determining the onset of convection that cover the likely range 
of convective core sizes for non-rotating stellar models — from 
Ledoux criterion to Schwarzschild criterion with convective 
overshoot. The qualitative conclusions of our stability analyses are 
insensitive to our choice of mixing parameters. Thus we present, as 
an illustrative example, models using the Schwarzschild criterion 
for convection and with convective overshoot of 10% of one 
pressure scale height above the c onvective core (with overshoot 
calculated using the prescription of lHerwi3l2000l) . 

We follow the evolution of Population III stars with masses 
between 120 and 1000 Mq through their main sequence (core hy- 
drogen burning) evolution, using the parameters outlined above and 
a nuclear reaction network coverin g all the relevant r eactions for 
hydrogen and helium burning (see Paxton et al.ll201 ll for details). 
Figure[T]shows the evolution of these stellar models from the Zero- 
Age Main Sequence (ZAMS) to a 40% central hydrogen mass frac- 
tion {Xcr). For massive Population III stars (M > 20 Mq), hydrogen 
burning by the p-p chain cannot halt the gravitational contraction; 
the star contracts to ~ 10** K and ignites He fusion by the triple- 
alpha process to produce enough heavy elements to s ustain hydro- 
gen burning by the CNO cycle jMarigo et aDl200ll) . Thus, these 
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Figure 1. Hertzsprang-Russell Diagram for Population III Stars with 
masses between 120 and 1000 and initially primordial abundances. For 
each mass, the evolution is run with Schwarzschild mixing and 10% con- 
vective overshoot, no rotation, and zero mass loss. Each sequence shown 
runs from the ZAMS to a center hydrogen composition (Xcn ) of 40% by 
mass. Filled symbols indicate models with unstable radial pulsations when 
convective damping is ignored, while open symbols denote stable models. 



stars have much higher central temperatures on the main sequence 
than their higher metallicity counteiparts. 

The highly mixed models shown in Fig.[T]have large convec- 
tive cores and correspondingly thin radiative envelopes, which ex- 
pand to form red supergiants while the core is still fusing hydrogen. 
This is particularly true for the most massive models, as can be seen 
in Fig. [10 



3 PULSATIONAL ANALYSIS 

We assume the quasi-adiabatic approximation for our pulsational 
analysis. For each model in the sequences shown in Fig. [T] (and 
the analogous higher metallicity models), we have calculated a set 
of adiabatic oscillation modes usin g the recently released version 
of the ADIPLS oscillation package JChristensen-Dalsgaardl private 



Evolving the most massive models to a lower central hydrogen mass 
fraction becomes difficult because the envelope expands significantly and 
evolves on a timescale much shorter than the nuclear timescale. For the less 
mixed models we have ran (those without overshoot or using the Ledoux 
criterion for convection), the radiative envelopes are larger for a given mass, 
and they do not expand to form red supergiants while still burning hydro- 
gen in their cores. Our qualitative results regarding the stability of main 
sequence massive stars are insensitive to our choice of mixing parameters. 
However, exactly when (and if) the star evolves significantly away from its 
ZAMS state does depend on the details of mixing. 
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communication). We assume zero Lagrangian pressure perturba- 
tion, 6p, as our outer boundary condition for all mode calculations. 
In all our calculations below, we use Sx to represent the Lagrangian 
perturbation of variable x, and x' its Eulerian perturbation. 

We calculate the perturbative driving an d damping rates fo r 
the adiabatic modes following the formalism of lUnno et al.l ( ll98^ . 
We thus approximate the imaginary component of the mode fre- 
quency (driving or damping rate) from the second-order expression 
for the entropy perturbation 



- ^ |(5rp dm 



(1) 



The above represents the integration of the change in entropy due 
to the perturbation (mode) over one oscillation period, 2n/a), and 
over the adiabatic mass, M„^i, of the star, defined as where the ther- 
mal time of the overlying mass is longer than the mode period (cf., 
[Noels 1998,, and our eq.ll4t. In this formalism, positive (negative) 
imaginary values represent mode driving (damping). Here 

^net con- 
tains the energy generation due to nuclear processes, as well as 
cooling from neutrinos (negligible here) and gravitational power 
from contraction or expansion (negligible here). 

We furth er separate the int egral in the numerator into work 
integrals, as in lUnno et al.l ( [T989h , such that 
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Equation|6]follows from* Aizenman & Coxl ( Il97 5h. but never makes 
any significant contribution for the modes investigated here since 

Cnuc on the main sequence. Note that we also ignore the effect 
of neutrino losses to the damping of modes since, again, fneut <K fnuc 
on the main sequence. 

Below we provide more details on the meanings of each of 
these work integrals. 

3.1 Nuclear Driving 

Oscillation modes with large amplitudes near regions of nuclear 
burning in the stellar interior can couple to this energy generation 
for mode driving. Over one pulsation cycle, ignoring other effects, a 
mode will experience a net gain in energy due to the excess energy 
generated during the phase of positive temperature perturbation, 
due to the strong, positive temperature dependence of most nuclear 
processes, d In e/d In T\p ~ 10 for hydrogen burning. 
Taking equation^ and plugging in for Se we find 
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where the sum over consecutive indices, /, represents a sum over the 
nuclear reactions taking place. The factors sj- and 6p represent the 



logarithmic partial derivatives of the energy generation rate against 
T and p respectively. We have also used the fact that the Lagrangian 
perturbations are adiabatic to eliminate the density perturbation in 

favor of the temperatur e perturbatio n. 

As pointed out bv lUnno et al.l jl989l) the values of Er,p are in 
general, frequency dependent. When the timescales associated with 
the individual reactions are sufficiently short or long compared to 
the mode period, simple approximations can be made to give appro- 
priate, modified € tm values (as in the treatment of the p-p chain in 
lUnno et alj[l989l) . However, when the timescales are comparable, 
one must account for phase sh ifts between the thermodynamic and 
abundance perturbations (cf., iKawalejl 19881 : ISonoi & ShibahashU 
[20I2V 

For the high central temperatures reached in Population III 
stars (logTr > 8.1), the fundamental mode period is inteiTnediate 
between the proton capture and beta-decay timescales in the CNO 
cycle. Accounting for the phases of the abundance variations, we 
find a slight reduction of fj- by ~ 15% from the equilibrium value 
determined by the temperature dependence of the '''N(/:), 7)'^0 re- 
action. Thus we have er ^ 8 near the centers of massive Po pulatio n 
III sta rs, somewhat smaller than the value used in iBarafFe et all 's 
( I2OOII) earlier study. 



3.2 Radiative Damping 

Equation|5]can be rewritten as 



Jo 



6T* d(SU) 
T dr 
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(9) 



where Lj^ signifies the radiative luminosity, and R„j is the radius 
corresponding to the adiabatic mass, M^rf, introduced in equation 
[T] above. Note that, according to this expression, damping occurs 
whenever the temperature perturbation and radial derivative of the 
luminosity perturbation are of the same sign. That is to say, if, at 
the temperature maximum, there is more perturbed luminosity leav- 
ing the top of the layer than entering the bottom, there is radiative 
damping. 

The luminosity perturbation itself has several terms; keeping 
only those relevant for radial oscillations, it is given by 



Lr r T k 
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where cj^ is the stellar natural frequency GMf,IR\, q is the local 
mass coordinate, m,-IMt,\ x is r/i?*; is the amplitude of radial 
displacement; V* is the logarithmic temperature derivative in the 
background model, dXnTldlnP; Wad is the adiabatic derivative, 
(31n r/31nP)s ; k the opacity; p the density; g the local gravity; 
and O' the Eulerian perturbation to the stellar potential. 

For oscillation modes with significant amplitude in the stel- 
lar core, radiative diffusion generally leads to modest damping 
of modes. Radiative damping becomes increasingly important as 
mode energy shifts into the stellar envelope. At low enough effec- 
tive temperatures and high enough metallicities, the opacity per- 
turbation term (Sk/k) can provide strong driving through the k- 
mechanism by blocking flux at compression (maximum positive 
temperature perturbation) such that d(SLR)/dr changes sign (e.g., 
.Cox.l98a) . 
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3.3 Convective Damping 

The characteristic convective damping ti mescale for stars with con- 
vective cores is ~ (M-i,R\lL-i,yi^ (e.g.. IZahrilll989l) . For massive 
stars, this is ~ 4 orders of magnitude shorter than the thermal time, 
which is the characteristic growth time of perturbations due to the 
e-mechanism (/c-driving can give rise to much more rapid growth). 
This highlights the importance of including convective damping in 
calculations of the pulsational stability of massive stars. The sub- 
tlety is that the effective convective viscosity is significantly sup- 
pressed relative to the above estimate because of the mismatch be- 
tween the convective turn over time and the mode periods of interest 
jGoldreich & Nicholsonl fig??: Zahn 1989). 

Recent simulations by Penev et al. (2009 a|) have s hown that, 
in agreemen t with the work of Goldreich & NicholsorJ (Il977h and 
IZahnHl989l) . interaction between pulsations and convective eddies 
indeed acts like a viscosit y, damping o scillatory modes in convec- 
tion zones. In particular, iPenev et al. I (|^09a) find that, over the 
range of frequencies and spatial scales accessible in their simu- 
lations, the effective kinematic viscosity of the convection scales 
linearly with the ratio of the oscillation period and turnover time of 
the largest eddies. We refer to the latter as Tl below. 

However, their simulations do not allow resolution of eddies 
with turnover times smaller than about 0.3 Tl- That they see no 
measurable viscous damping for forcing periods less than this res- 
olution limit implies that eddies with turnover times less than or 
equal to the forcing timescale dominate the viscous interaction 
w hen fforce i 0.3 T/^. Thus, these authors argue that the formalism 
of iGoldreich & Nicholsoiil ( Il977l) . derived on the assumption that 
near-resonant eddies in a Kolmogorov cascade dominate the vis- 
cosity, applies for higher frequency (shorter period) forcing. 

We thus take a conservative approach, in terms of minimiz- 
ing the effect of viscous damping, and assume that the change to 
quadratic scaling occurs just beyond the resolution of their simula- 
tio ns, at a ratio of far c e /tl = T^mm = 0. 1 . Generalizing the results 
of IPenev & Sasselov ( |20T l|), we define our fiducial, frequency- 
dependent kinematic viscosity to be 



v(<jj,r) = I -VLL|min 



1 



2n 



In 



(11) 



Here r is the radius of the layer, L = min(r, //) is the mixing 
length, vl is the convective velocity according to mixing length the- 
ory, Ti is the convective turnover time of the largest eddies (which 
we take to be 1 /A^Bmnt-vaisaia), w is the mode's angular frequency 
(and 27r/a) = /force), and the dimensionless factor H„,ax = 2.4 
represents the timescale ratio above which the viscosity is inde- 
pendent of forcing period (saturation period). As can be seen, the 
quadratic reduction applies for (2;r/wT/.) < Hmin, linear reduction 
for H„,i„ ^ {InlujTi) < H„,„, and unreduced for (2n/cjTi) > H„,„ 

Figure|2]shows a comparison of Tl and the fundamental mode 
period for several main sequence models. The left panel shows Pop- 
ulation III models, while the right shows higher metallicity, 120 Mq 
models. For Population III models, the fundamental mode periods 
are always well into the quadratic scaling regime, but for the higher 
metallicity cases, the mode periods tend to lie in the linear scaling 
or unreduced regimes, especially for the convective zones in the 
envelope of the star (log T < 6). 

For a radial mode within the star, the integrated rate of con- 
vective damping is given by 



dm. 



(12) 
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Figure 3. Growth rates, in units of Myr" ' , for the radial fundamental modes 
of Population III stars on the main sequence, shown versus the mass fraction 
of hydrogen in the star's convective core. The top panel shows the results as- 
suming the "frozen-in" approximation for convection (for comparison with 
fig. 4 of Baraffe et al. 2001), while the bottom panel includes the effects of 
convective damping. Negative growth rates (shaded region) denote stabil- 
ity. Note that the 120 M© models are always stable in this analysis with 
or without convective damping, due to the lower CNO cycle temperature 
dependence (see Section lxTt 



Plugging in the expansion in terms of modes, we find 



£'co„v = J Vs'qIiJ^ \k,f ^] 



dm. 



(13) 



where s'^ is the coefficient for damping due to the radial component 
of the shear from Penev & Sasselov (2011, their eq. 9), and k, is 
given by the usual adiabatic oscillation equations. In accordance 
with equation [T] we then divide by the mode energy to calculate 
the integrated damping rate due to interaction with the convection. 

For the modes in Population III stars described below, 
(Inlwrt) < Tlmi„. If the turnover to quadratic scaling actually oc- 
curs at mode periods shorter (longer) than what we have taken 
here, the damping will be stronger (weaker) than what we calcu- 
late. However, give n the limiting value of n,„i„ < 0.3 found by 
IPenev et "ah (2009a). we find it unlikely that the damping is weaker 
than calculated here (see Sections|4]and|5]for further discussion). 



4 RESULTS 

4.1 Zero Metallicity 

Figure [3] shows the results of our quasi-adiabatic stability analysis 
for the fundamental radial modes in Population III main sequence 
models. The top panel shows the driving/damping rates using the 
"frozen-in" approximation for convection, while the bottom panel 
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Figure 2. Convective turnover time, t/,, plotted against fractional mass enclosed for a Population III main sequence model with masses from 120 - 10^^ Mq 
(left panel), and against internal temperature for 120 M© models for a range of metallicities (right panel). For the Population III stars, the models shown have 
Xctr ~ 0.7; the vaiiable metalhcity stars have X^r ~ 0.4. The straight horizontal lines in each panel show the fundamental mode periods coiTesponding to each 
model plotted. For reference, the boundary between linear and quadratic scaling of the convective viscosity with mode period for the IOOOMq Population III 
model is shown as the thick dashed grey line in the left panel. For all other models (in both panels) this transition occurs at a period one order of magnitude 
less than the convective turnover time of the stellar model. 



includes the effects of convective damping. In the absence of con- 
vective damping, massive stellar models with M > 120 Mq have 
unstable fundamental modes until the central hydrogen mass frac- 
tion is < 0.6. The bottom panel of Fig. [3] demonstrates, however, 
that convective damping stabilizes the fundamental mode in all 
Population III main sequence models, regardless of mass and age. 
We find that this is true unless n„,i„ (which characterizes the tran- 
sition from linear to quadratic suppression of convective damping) 
satisfies Tlmm <; 0.25. While this possibility is not completely ruled 
out by Penev et al.'s calibration of convective damping, we regard 
it as exceptionally fine-tuned. 

When we ignore the effects of co nvective damping (to p panel), 
our results are comparable to those of lBaraffe etalj toll their fig. 
4) to within factor of 2. The moderate disagreement likely stems 
from differences in the parameters used to evolve the stellar mod- 
els and the ~ 15% difference in the temperature dependence of the 
CNO cycle between their work and ours (see Section l3Tt . In par- 
ticular, the reduction in temperature sensitivity leads us to conclude 
that the fundamental mode in 120 Mq Population III main sequence 
models is stabilized by radiative diffusion in the envelope even be- 
fore convective viscosity is taken into account. 

While our calculations are quasi-adiabatic and theirs are non- 
adiabatic, this difference is likely unimportant for the modes and 
models shown in Fig. [3] which Baraffe et al. found to be unstable to 
the 6-mechanism (stars having > 0.5). First, the non-adiabatic 
region of the star, defined as where 

t,, ['^ Cv T dm 

'""="" = < 1 (14) 

^mode ^mode 

contains less than ~ lO""* of the total stellar mass for all models 
shown in Fig. [3] Second, the fraction of the mode energy found 
in the non-adiabatic zone is also always less than ~ lO""*. Taken 
together, these justify our quasi-adiabatic assumption. 

Figure |4] shows the mode energy distribution (top panel) and 
work integrals (bottom panel) for the fundamental radial mode in 
a main sequence (JC^,. = 0.718) 500 Mq, Population III stellar 
model. The dominant contributions to the total work integral come 
from the nuclear driving and convective damping, at r/R < 0.4 
(fMr :£ 100 Mq). Similar plots for fundamental modes in other main 
sequence models are qualitatively the same. For more evolved mod- 
els (those with lower Xar), the mode energy shifts outward into the 
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Figure 4. Mode energy (red solid line, top) and work integrals (bottom) 
for a typical, damped radial fundamental mode in a Population 111, SOOMq 
main sequence model. All quantities are plotted versus r/R (bottom abcissa) 
and enclosed mass, (top abcissa). The top panel also shows the nuclear 
energy generation rate (blue dashed, right axis) and the Mach number of 
the convection (purple, dotted line). The bottom panel shows the different 
contributions to the total work integral described in equations |4] ID [H and 
1131 with colors according to the legend. 

stellar envelope and radiative diffusion is the dominant source of 
damping. 

4.2 Dependence on Metallicity 

We have performed the same analysis as for the Population III stars, 
for models with 120 Mq and a range of metallicities from 10"* to 
solar (taken to be 0.02). Figure |5] shows the results of these calcu- 
lations. In the upper panel, the net growth rate (taking into account 
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Figure 5. Net growth rates including convective damping (top) and the con- 
vective damping rates (bottom) for radial fundamental modes in stars on the 
main sequence with M = 120 Mq and a range of metallicities, shown ver- 
sus the mass fraction of hydrogen in the star's convective core. Note the 
logarithmic scale of the y-axis, which highlights the large change in the 
magnitude of the driving/damping when the mode amplitude shifts to the 
envelope, as for the solar and half solar metallicity models (blue circles and 
green squares, respectively). The corresponding change in the magnitude 
of the convective damping occurs because the damping is dominated by an 
envelope convection zone rather than the convective core (see Fig.|2). 



ment changes with the expansion of the stellar envelope, such that 
only models at particular evolutionary stages on the main sequence 
have fundamental modes driven by the /c-mechanism. 

Across the range of masses and metallicities investigated here, 
we find that convective viscosity always overcomes the e-driving in 
the stellar core. However, for models with Z > 2 x 10"', the stel- 
lar and mode structures are such that /c-mechanism provides the 
dominant driving and is orders of magnitude stronger than the con- 
vective damping. However, this envelope driving occurs in a re- 
gion where the local ratio f therm /f mode is only slightly larger than 
unity and our assumption of quasi-adiabaticity is not entirely valid 
(as shown in Fig. |6j. A fully-non-adiabatic method would likely 
provide more reliable growth-rates for these modes. Moreover, our 
quasi-adiabatic analysis is not capable of capturing strange-mode 
driving, which is also likely to be impor tant for the sa me models 
that are unstable to the /c-mechanism (cf.. lGlatzeill 1998h . 

4.3 Higher order and non-radial modes 

In addition to the fundamental mode, we have also studied the 
stability of higher order radial modes and non-radial oscillations. 
Higher order radial modes — those having a larger number of radial 
nodes — have a more significant fraction of their energy outside the 
nuclear burning core, and have shorter wavelengths, which both 
lead to stabilization by radiative damping. 

The presence of a large convective core prevents gravity (g)- 
modes from pro pagating to small radii and coupling to the nuclear 
driving (see also lSonoi & Umedal201 ll) . While high frequency, low 
angular degree pressure(p)-modes may penetrate into the burning 
region, these modes have their largest amplitudes near the stellar 
surface where the density is low and radiative damping dominates. 
Furthermore, while p- and g-modes of high enough angular degree 
can be trapped by the increase in A'^ at the edge of the convective 
core due to the chemical composition gradient, these modes expe- 
rience strong radiative dam ping that overcom es the weak nuclear- 
driving they experience (cf. lUnno et al 



any convective damping) is shown versus the central hydrogen frac- 
tion (as in Fig. [3]l. Note the logarithmic scale for the growth and 
damping rates. The bottom panel shows the contribution of the con- 
vective viscosity to the damping rate. Models with < Z < 2x lO""* 
are not shown since their results are nearly identical to the Z = 
models. 

As shown in Fig. [6] the sharp transition in the net growth and 
convective damping rates for models with Z > 2 x 10"' is due to 
the shift in mode energy out into the stellar envelope (shown in the 
figure for models with Z = 0.01). This transition occurs when the 
low-density envelope expands due to the increased radiation pres- 
sure around the iron opacity enhancement at log T x 5.2. When the 
mode energy is primarily in the stellar envelope there is the pos- 
sibility for /c-mechanism driving around the iron opacity enhance- 
ment, as well as stronger convective damping in an outer convection 
zone where the mode a nd convective freque ncies are comparable. 

In agreement with lBarafi^e et alj ( l200lh . we find that the radial 
fundamental modes in main sequence models with Z > 2 x 10"' 
are strongly driven by the /c-mechanism operating in the envelope 
at the locat ion of the opacity-enhan cement due (primarily) to iron 
transitions jRogers & Iglesiaslll992[) . The location of the adiabatic 
cutoff mass/radius, as defined by equation [14] relative to the loca- 
tion of the opacity peak at log T x 5.2 determines when the opacity 
is capable of providing strong driving (cf. iCoxiil980i) . This align- 



5 DISCUSSION AND CONCLUSIONS 

In this paper, we have studied the stability of radial and non-radial 
oscillation modes for massive main sequence stars of a range of 
metallicities, from primordial Population III stars to those of solar 
composition. We find that convective damping is likely to stabi- 
lize modes that would otherwise be unstable by the e-mechanism 
(driving by nuclear fusion). This is particularly important for pri- 
mordial stars whose high effective temperatures imply that there is 
no driving due to opacity variations in the stellar envelope (the k- 
mechanism). Thus, we conclude that massive Population III stars 
are linearly stable on the main sequence, supporting the idea that 
such stars are unlikely to experience any pulsation-driven mass loss 

during their main sequence evolution. 

The c alculations in this paper build on those of ISarafi^e et al.l 

l l200lh and lSonoi & Umedal ( 1201 ih , who found that the fundamen- 
tal radial oscillations in Population III stars with M > 120 Mm are 
unstable due to the e-mechanism. However, both lSarafi^e et al.l and 
Sonoi & Umeda neg lected convective da mping. Recent hydrody- 
namic simulations bv lPenev et al. U2009ah have shown that the con- 
vective viscosity scales linearly with the oscillation period for peri- 
ods of the same order as the turnover time of the largest eddies (in 
agreement with Zahn 1989) and quadratically for much shorter pe- 
riods (in agreement with ,Goldreich & Nicholson 1977). The tran- 
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log T [K] log T [K] 

Figure 6. Mode energy (red) as a function of log T for Z = 0.01 main sequence stellar models. Shown on the left is a stable mode with Xdr ~ 0.58, and on 
the right is a Ar-driven mode with X^r ~ 0.48. Also shown are the run of opacity (green dot-dashed line, left ordinate, multiplied by 3 for clarity), convective 
mach number (purple dotted line, right ordinate), nuclear energy generation rate (blue dashed line, left ordinate), and adiabatic cutoff (black dashed vertical 
line). All values are shown versus log T. In the more evolved model on the right, the low-density envelope has expanded by a factor of ~ 2 due to the enhanced 
radiation pressure around the iron opacity bump at log T » 5.2. 



sition between these regimes is not directly accessible to the sim- 
ulations; we have taken a conservative estimate of n,ni„ = 0.1, in 
agreement with the simulations, in order to minimize the convective 
viscosity. If ITniin > 0.25, some e-driven mo des remain u nstable, but 
this is an extremely fine-tuned value, given lPenev et alj s result that 

Hmin < 0.3. 

Even more recent simulation work by lOgilvie & Lesuj 

( l2012h calls the work of IPenev et^ into question, finding that 
the viscosity always scales with the period ratio squared (as 
iGoldreich & Nicholson|[T977l) . and may in fact be negative (i.e., 
driving) in the short period forcing regime of interest for this work. 
T hese latter sim ulations use a smaller forcing amplitude than that 
of IPenev et al.L leading to relatively noisier measurements, and a 
smaller simulation box which may not capture the full spectrum 
of turbulence on the large scales within a stellar convection zone. 
This latter aspect may be crucial to their measurement since they 
infer that the negative contribution to the viscosity is dominated by 
the effect of the largest scale eddies. The sensitivity of our results 
to the exact form of the convective viscosity highlights the need 
for higher resolution simulations of the convection-oscillation in- 
teraction, particularly using larger simulation domains to capture 
the interaction with the large-scale eddies. 

Stellar models with sufficiently high metallicity, Z > 2 x 10"', 
have strong opacity variations in their envelopes which can provide 
/(-mechanism driving. We find that the fundamental mode in these 
models are destabilized by the ^c-mechanism operating near the 
iron-bump in opacity for much of their main sequence evolution. 
This driving is orders of magnitude stronger than the convective 
damping. These models experience large linear growth rates that 
approach non-adiabiticity (5a) < to). For lower metallicity models, 
the e-mechanism provides the dominant driving but is overcome 
by convective damping as in Population III models. In all of the 
massive stellar models we have considered, we find that damping 
due to convective viscosity is stronger than the driving by the 6- 
mechanism in stellar interiors for radial modes. 

Both lBaraffe et alj and lSonoi & Umed3 estimate that massive 
Population III stars may lose a few (for M < 1000 Mq) or up 
to ~ 10% (for M > 1000 Mg) of their total stellar mass due to 
pulsation-driven mass l oss during the pr oposed unstable phase of 
the main sequence. As ISonoi & Umedal highlight, the convective 



cores of stars > 500 Mq comprise > 90% of the stellar mass, imply- 
ing that pulsation-driven mass loss may expose the convective core, 
perhaps further increasing mass loss by mixing nuclear-processed 
material to the surface to initiate a radiation-driven wind, or by in- 
creasing the duration of the proposed pulsationally unstable phase. 
If, instead, the pulsations are damped by convective viscosity as we 
have concluded, there may be no significant mass lost during the 
main sequence phase for massive Population III stars. 

In order to fully address the mass loss properties of massive 
primordial stars, we must address their stability at more evolved 
stages. However, as stars evolve off the main sequence, their fun- 
damental radial mode becomes increasingly non-adiabatic as the 
mode energy shifts to the expanding, low-density stellar envelope. 
Analyzing these modes requires a non-adiabatic calculation, which 
we leave to future work. 

With all oscillatory modes in main sequence Population III 
stars likely stabilized by convective damping, evidence points to 
there being no mass loss on the main sequence for these metal- 
free stars, so long as the non-rotating stellar models used here are a 
reasonable approximation. The question of what happens as these 
stars evolve thus becomes all the more crucial for assessing the 
fates of the first stars. 
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